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We consider the cosmological consequences if a small fraction (/ < 0.1) of the dark 
matter is ultra-strongly self-interacting, with an elastic self-interaction cross-section 
per unit mass cr 3> 1 cm^/g. This possibility evades all current constraints that 
assume that the self-interacting component makes up the majority of the dark mat¬ 
ter. Nevertheless, even a small fraction of ultra-strongly self-interacting dark matter 
(uSIDM) can have observable consequences on astrophysical scales. In particular, 
the uSIDM subcomponent can undergo gravothermal collapse and form seed black 
holes in the center of a halo. These seed black holes, which form within several 
hundred halo interaction times, contain a few percent of the total uSIDM mass in 
the halo. For reasonable values of af, these black holes can form at high enough 
redshifts to grow to ~ 10® Mq quasars hy z > 6, alleviating tension within the stan¬ 
dard ACDM cosmology. The ubiquitous formation of central black holes in halos 
could also create cores in dwarf galaxies by ejecting matter during binary black hole 
mergers, potentially resolving the “too big to fail” problem. 
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1. INTRODUCTION 

Although ACDM cosmology provides an excellent fit to the observational data on > Mpc 
scales [1], its success is less certain over the strongly nonlinear, < kpc regime relevant to the 
substructure within galactic halos. The deviation of galactic cores from an expected cuspy 
density profile Eli and an apparent shortfall of observed Milky Way satellites relative to 
expectations from simulations Hi originally motivated considerations that the dark matter 
might have non-negligible self-interactions [B]. Although a combination of improved theo¬ 
retical understanding and additional observations had appeared to alleviate these problems 
and remove the phenomenologically interesting parameter space for self-interacting dark 
matter (SIDM) [3-[9], a recent reevaluation of the constraints m ITT] has demonstrated 
that SIDM with an velocity-independent elastic self-interaction cross section per unit mass 
cr ~ 0.1 — 1 cm^/g ~ 0.2 — 2 b/GeV can simultaneously meet all constraints and alleviate 
the discrepancies between ACDM and observations. 

In this paper, we exhibit a distinct area of the SIDM parameter space which is likewise 
both allowed by observations and potentially interesting phenomenologically. In particular, 
we examine the case in which most of the dark matter remains non-self-interacting (or 
weakly self-interacting) as in the standard ACDM picture, but a small fraction / -C 1 of the 
dark matter is made up of a subdominant component that is ultra-strongly self-interacting, 
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abbreviated as uSIDM, with a 3> 1 cm^/g (where a denotes the cross section per unit mass). 
Because most of the dark matter remains inert, constraints that rely on distinguishing the 
overall behavior of SIDM halos from their CDM counterparts are no longer relevant. 

Consider, for example, the constraints placed on the SIDM cross section from observations 
of the Bullet Cluster (IE 0657-6). Observations reveal an offset between the gas “bullet” and 
the dark matter centroid of the currently merging subcluster. Under the assumption that 
the subcluster has already passed through the main cluster, this offset is due to stripping 
and deceleration of gas in the subcluster due to interactions with the main cluster itself. 
The observation that the dark matter has not been slowed to the same degree allows limits 
to be placed on the dark matter self-interaction cross section. The strongest constraint [12] 
comes from the measurement of the ratio of mass-to-light ratios of the subcluster and the 
main cluster, which is found to be 0.84 ± 0.07. Under the assumption that the subcluster 
and main cluster had the same initial mass-to-light ratio before merger, this means that the 
subcluster cannot have lost more than 23% of its mass. 

In [I2|, this measurement plus estimates of the subcluster escape velocity and merger 
speed were used to constrain a < 0.6 cm^/g when / = 1. However, it is clear that, even 
in the extreme example that all of the SIDM mass in the subcluster was lost to scattering, 
current observations would not be able to detect the SIDM subcomponent if / < 0.07, within 
the uncertainty on the mass-to-light ratio. So constraints from the Bullet Cluster certainly 
do not apply when / -C 10“^, regardless of the size of the self-interaction cross section per 
unit mass a. Even when / ~ 0.1, a may not be well-constrained, since one pass through 
the main cluster would not suffice to strip all of the SIDM from the bullet. 

We note that observations of another cluster undergoing a major merger, A520 unHini 
have not provided similar constraints on the SIDM cross section; here the dark matter 
centroid of the subcluster is in fact coincident with the (presumably stripped) gas. Under 
certain assumptions, this can be taken as evidence of a nonzero dark matter self-interaction 
cross section per unit mass, as strong as 0.94 ± 0.06 cm^/g in the latest observations [T6] . 
The limited number of ongoing major merger events in the observable universe makes it 
hard to give an overall estimate of the self-interaction cross section from major mergers, but 
future surveys could potentially combine many minor merger events to measure a with a 
precision of 0.1 cm^/g [T7] . 

Regardless of the situation for / = 1 SIDM, we have seen that there are no observational 
constraints on a uSDIM component of the dark matter with / < 0.1. At the same time, 
of course, a small component of uSIDM by itself is unable to produce cores or dissolve 
substructure to any observable degree. We point out, though, that a uSIDM component 
of the dark matter could instead explain another potential discrepancy with the ACDM 
picture: the existence of billion-solar-mass quasars at high redshifts z > 6.5 — 7 (for reviews, 
see [THH2SI)- In §2 we review the observational situation and the difficulties with explaining 
it within ACDM. In §3 we suggest an alternative: gravothermal collapse of an ultra-strongly 
self-interacting dark matter component. We review the mechanism of gravothermal collapse, 
specialize to the case of a halo containing uSIDM, and solve the problem numerically. We 
apply the results of §3 to individual observations of high-redshift quasars in §4, then discuss 
broader cosmological implications in §5, including a potential way for uSIDM to indirectly 
produce cores in dwarf halos. We hnally conclude in §6. 
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2. SUPERMASSIVE BLACK HOLES 


Supermassive black holes (SMBHs) which grow primarily via gas accretion are Eddington- 
limited: the gravitational force on the accreting gas is balanced by its own radiation pressnre. 
Hence growth via gas accretion cannot proceed faster than exponentially, with an e-folding 
rate bonnded by the inverse of the Salpeter time m- 
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ATTGrrir, 


^ ) 45.1 Myr, 


where ctt is the Thompson cross section, 
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rup and me are respectively the proton and electron masses, and e^ is the radiative efficiency, 
which ranges from 1 — -^8/9 ~ 0.057 to 1 — -^1/3 ~ 0.42 as the angnlar momentum of the 
black hole increases from zero to its extremal value [2S]; in astrophysical applications, e,, is 
typically taken to be = 0.1. Accretion faster than the Eddington limit, MEdd = ^^sai) 
onto a black hole of mass M will result in a radiation pressure exceeding the gravitational 
force, driving outflows which should quickly halt this excessive accretion. Yet several dozen 
quasars with masses a few x 10® M© have been detected at redshifts ^ > 6, including 
a quasar, ULAS J1120-I-0641, with mass 2.0^1'^ x 10® M© at redshift x = 7.085 [26l |27]. 
Using the Planck Collaboration’s best-£t cosmological values [I], x = 7.085 corresponds to 
747 Myr after the Big Bang, so even continuous Eddington accretion since the Big Bang 
can only increase the mass of a seed black hole by a factor of 1.6 x 10^. If we make the 
standard assumption that black hole seeds are formed from Pop III stars, the seed cannot 
have formed before around x ~ 30, so the maximum growth factor shrinks by another order 
of magnitude, to 1.75 x 10®, requiring a seed black hole mass ~ lO^M©. 

More generally, in order to explain the observed abundance of ~ 1/Gpc^ billion-solar- 
mass quasars at x ~ 6 [23] within ACDM, we must form 10^“^M© seed black holes soon 
after the beginning of baryonic structure formation and grow these black holes continuously 
at near-Eddington rates for ~ 800 Myr. Some simulations have shown this can be achieved 
[2H|, but only by making optimistic assumptions about cooling and star formation [221 ED], 
fragmentation [3TVI33] . photoevacuation [34113B] . black hole spin [571139] . and black hole 
mergers [401 - 142] . We emphasize, in particular, that these results depend critically on the 
assumption of = 0.1; because the e-folding time itself depends linearly on the radiative 
efficiency, the maximum mass formed by a given time is exponentially sensitive to its value. 
Because quasar masses are inferred by measuring their luminosities and assuming they are 
Eddington-limited, increasing the assumed radiative efficiency will decrease the inferred 
quasar mass by e~^. However, this reduction in required mass is made negligible by the 
much larger number of e-folds required to reach it. Recent work, both theoretical [22] and 
observational ]43], has found ^ 0.2, which would be catastrophically incompatible with 
an assumption of black hole growth driven by Eddington accretion. 

One alternative is to allow for extended periods of super-Eddington gas accretion. Super- 
Eddington accretion is known to be possible, for example when outflows of gas and radiation 
are collimated [441145] . and extended periods of super-Eddington growth could account for 
the observed supermassive high-redshift quasars [ig HZ]. However, estimates of quasar 
masses and luminosities at low redshifts using emission line widths indicate that, at least in 
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the late universe, the vast majority of quasars are constrained to radiate at the Eddington 
limit [H], or possibly well below it |l9l 150] . 

In this paper, we will therefore neglect the possibility of extended super-Eddington ac¬ 
cretion. We will assume that growth of black holes from baryonic accretion is limited to 
exponential growth with an e-folding time given by the Salpeter time 0 . In order to facil¬ 
itate comparision of uSIDM to the standard picture, we will, however, allow for continuous 
accretion of baryons at this limit once a seed black hole has formed, despite the potential 
issues mentioned in the previous paragraph. In other words, we attempt to modify the 
mechanism by which black hole seeds are formed, while leaving the simplest conventional 
mechanism for their growth from seeds to supermassive black holes intact. It would be easy 
to combine our results with more realistic baryon accretion histories. 

Finally, we note that future observations in the near-infrared, e.g. with the James Webb 
Space Telescope and Wide Field Infrared Survey Telescope (WFIRST), and in the radio, e.g. 
with the Square Kilometer Array, should be able to detect (or place limits on the density 
of) even intermediate-mass (~ lO^M©) quasars out to z ~ 10 [5TH51], providing vastly more 
information about the formation and growth of high-redshift quasars. 


3. GRAVOTHERMAL COLLAPSE 

Motivated by the tensions within the standard (ACDM) picture discussed in the previous 
section, we propose an alternative mechanism for black hole seed formation: the gravother- 
mal collapse of the uSIDM component of a dark matter galactic halo. The simplest 
form of gravothermal collapse occurs in a population of gravitating point particles with elas¬ 
tic short-range interactions. The classic illustration of the mechanism is globular clusters, 
where the point particles are stars. Stellar short-range interactions are not purely elastic, 
so in this case collapse is eventually halted by binary formation. A gas of SIDM, however, 
has only elastic interactions, so core collapse continues until relativistic instability results in 
the formation of a black hole, which promptly Bondi accretes |56] the optically thick core 
of SIDM that surrounds it. 

In this section we make this intuitive picture precise. Full expressions will be given below, 
but in brief we hnd that the uSIDM component of a galactic halo undergoes gravothermal 
collapse in ~ 460 halo relaxation times, forming a black hole which contains ~ 2% of the 
uSIDM mass of the galaxy. The halo relaxation time is a complicated expression which 
depends on the halo mass and time of formation as well as the uSIDM properties, but we 
show in the following sections that, for reasonable values of uSIDM fraction / and cross 
section per unit mass a, there exist halos that can easily form seed black holes, and grow 
them using uSIDM and baryons, to achieve IO^Mq SMBHs by redshift 6. 

Before formulating the problem, we hrst review the gravothermal collapse mechanism 
itself. Intuitively, gravothermal collapse depends on the simple observation that gravitation¬ 
ally bound systems have negative specihc heat. For a virialized system, this is immediate: 

0 = 2T + V = T + E^E = -T. (3) 

Now consider two systems, an inner, gravitationally bound system with negative specihc 
heat and an outer system surrounding it with positive specihc heat—the inner and outer 
parts of a globular cluster, for example. Evolution towards equilibrium will direct both mass 
and heat outward, causing both the inner and the outer system to increase in temperature. 
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A possible physical mechanism is a two-body scattering in the inner system which sends one 
star closer to the core (where it gains potential energy and thus speeds up, increasing the 
temperature of the inner system) and kicks one star out to the periphery (where its higher 
speed increases the temperature of the outer system). Importantly, we see that the inner 
system shrinks as it heats up. 

Now two outcomes are possible, depending on the specihc heat of the two systems as a 
function of their masses. If the outer system always has the smaller (magnitude of) specihc 
heat, its temperature will eventually grow to exceed that of the inner system, and the entire 
assemblage of masses will reach equilibrium. On the other hand, if the outer system grows 
in mass too quickly, its specihc heat will become too large and its temperature will never 
catch up to the inner system. Hence the inner system will continue shrinking in mass and 
growing in temperature until the thermodynamic description breaks down. This is precisely 
the gravothermal catastrophe [55] . In the case of a globular cluster (at least an idealized one 
with uniform-mass stars), the gravothermal collapse process is halted by binary formation, 
which acts as an energy sink [57] [58] . If the uSIDM interacts purely via elastic scattering, 
however, no bound state formation is possible, and gravothermal collapse can drive the core 
to relativistic velocities, where it undergoes catastrophic collapse into a black hole via the 
radial instability [591 - 162] . 


A. The Gravothermal Fluid Equations 


We now consider the gravothermal collapse of a general two-component dark matter halo, 
where the self-interacting component comprises some fraction / of the mass of the halo. At 
this stage we do not yet specialize to the uSIDM case, with / <C 1. To avoid confusion, 
we will therefore refer to the two different components of the halo as SIDM (making up a 
fraction / of the total mass of the halo) and (ordinary) CDM (making up the remainder), 
denoting the SIDM as uSIDM only when / <C 1. To simulate the collapse, we employ the 
gravothermal fluid approximation [631165] . which reduces the problem to a set of coupled 
partial differential equations that can then be solved numerically. First consider the general 
case for an / = 1 fluid, i.e. a halo composed entirely of SIDM. A spherically symmetric 
ideal gas of point particles in hydrostatic equilibrium with arbitrary conductivity k obeys 
the following equations [63] : 


dM . 

(4) 

—— = Airr^p 
or 

d (pz/2) GMp 

1^2 

(5) 

L dT 

(6) 

dr 

-dvrprV^In—, 

(7) 


where z/(r) is the one-dimensional velocity dispersion and L(r) the total heat radiated in¬ 
ward through a sphere of radius r. The hrst equation (|^ simply dehnes the integrated mass 
distribution in terms of the density. The second (|^ is the statement of hydrostatic equi¬ 
librium: we inserted Euler’s equation into the Poisson equation for a spherically symmetric 
potential and used the equation of state for an ideal gas, p = The third ([^ states that 
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the heat flux is proportional to the temperature gradient, with proportionality constant 
given by the conductivity k. The fourth ([^ is the second law of thermodynamics, inserting 
the specific entropy of an ideal gas of point particles u = ^ using the relation 

= k-QT/m. This gives a set of four differential equations with four dependent variables 
{M, p, z/, L} and two independent variables {r, t}. (The temperature T is directly related 
to u hy = k-QT/m.) 

To make progress, we need an expression for the form of the thermal conductivity k 
in terms of our physical parameter, the elastic scattering cross section per unit mass a. 
Dimensional analysis alone will not suffice: we have one time scale, the fluid relaxation time 


tr = l/(ap(TZ/), 


( 8 ) 


with a = ~ 2.257 for hard-sphere interactions, but two length scales, the mean 

free path A = l/{pa) and the Jeans length or gravitational scale height H = 

Following inn [65], we find the unique length scales in the two limiting cases, the short 
mean free path (smfp) regime A -C 77 —)■ 4mfp = and the long mean free path (Imfp) 
regime A 3 > 77 —)■ (.imfp = 77, and combine them in reciprocal to get a final length scale, 
7 = (7^fp -|- . In the smfp regime, transport theory tells us that 


7 3 

4^ “ “2“ 


(9) 


The coefficient b is calculated perturbatively in Chapman-Enskog theory | 66 ], b = 
25y^/32 cs 1.385. In the Imfp regime, the flux equation is well approximated as 


L _ _3 77^ 

A'lrr'^ 2 ^ dr ^ 


( 10 ) 


where C* is a constant setting the scale on which the two conduction mechanisms are equally 
effective, determined by N-body simulations [US] to be C ~ 290/385 ~ 0.75. Hence the final 
expression is 


7 

jTrr^ 
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b AttGI 
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Now consider the more general case, / 7 ^ 1. Hydrostatic equilibrium is separately satisfied 
for each species of particle, but the gravitational potential is of course sourced by both 
species, giving the coupling between the two components. Because the non-SIDM component 
is taken to be collisionless, it has a = 0, so 7™ = 0. So the total system is governed by 
six partial differential equations with six dependent variables {il7, p*"'*, p"'*, z/**^*, z/™, 7*"'^} 
and two independent variables {r, t}: 


^ (p‘"‘ + r) 

a (p“ (!/“)") 

dr 


( 12 ) 

(13) 


dr 


( 14 ) 
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(16) 

(17) 


As before, the first equation gives the total mass distribution, while the second and third 
enforce hydrostatic equilibrium. The fourth determines how the SIDM fluid conducts heat 
and the hfth how the flux gradient affects the fluid. Finally, the sixth equation ensures 
that the entropy of the collisionless component is conserved, 3h/z/ = p/p. Notice that the 
fraction / does not appear in the differential equations themselves, but only in the boundary 
conditions: we must have 


47rr'^p®"'*(r')(ir' / 

47rr'2p™(r')(ir' 1 — / 


(18) 


at all times. 


B. Initial Conditions 


In principle (12 -17) can be solved exactly given appropriate boundary conditions at 


r = 0 and r = oo and a set of initial radial prohles which obey the equations. In practice, 
this is computationally impossible: even hnding the initial prohles for an arbitrary a is 
infeasible. Balberg, Shapiro, and Inagaki [M], considering the / = 1 case, took the cr —)■ 0 
limit, which admits a self-similar solution where separation of variables is possible, then 
found the eigenvalues of the resulting system of ordinary spatial differential equations and 
took the resulting prohles as their initial conditions for the more general a ^ 0 case. 

We will instead assume that SIDM self-interactions are unimportant during the process 
of halo formation, so that the the SIDM and collisionless components have the same initial 
prohle. This allows us to use the results of (collisionless) ACDM simulations. We simplify 
further by approximating the initial halo by an NFW prohle. 


Pnfw(^ ) — 


Ps 


(r/rj(l -h r/rJ2’ 


(19) 


where ps and are the characteristic density and scale radius, respectively. Since the NFW 
prohle has a characteristic radius, we can state our assumption more precisely: we assume 
that halo formation proceeds much faster than heat conduction, which is true when the 
dynamical timescale of collapse is much less than the relaxation timescale due to collisions: 

tdynirs) < t^elirs) ~ “tdyn(?’s) -)■ <1; (20) 

A 

i.e. so long as the halo is optically thin at its characteristic radius. Again, if the optical 
depth is small, 

T = /Pnfw^^ct < 1 a/ < -, 

PsTs 


( 21 ) 
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typical SIDM particles have not yet undergone any self-interaction by the time of halo for¬ 
mation, so we are justified in assuming they follow the same initial prohle as the collisionless 
dark matter, = /pnfw(’")- 

Before checking the validity of this assumption, we comment on the consequences of taking 
a different initial prohle. The NFW prohle is particularly simple: its form means that the 
optical depth at small radii, r <C r^, is independent of radius, so a small characteristic 
optical depth implies that the central regions are also optically thin despite the presence 
of a cusp. Modern ACDM simulations, however, have tended to hnd density prohles more 
complicated than the NFW prohle. Prohles with cores or at least less cuspy behavior, e.g. 
Einasto prohles |63l68], will have r <C 1 everywhere if -C 1. Below we will see that SIDM 
halos with initial NFW density prohles grow cores on a scale of tens of halo relaxation 
times anyway, so shallower initial prohles will only result in slightly smaller times before 
black hole formation. Prohles with more cuspy behavior, e.g. generalized NFW or Zhao 
prohles [69j with inner slope a > 1, will unavoidably have regions at very small radii in the 
optically thick regime. Below we will see that SIDM halos with initial NFW prohles hrst 
evacuate the cusp to form cores before beginning the gravothermal collapse process, and it 
seems reasonable to conclude that the same thing will happen for non-pathological cuspier 
prohles. We conclude that imposing a diherent prohle should not signihcantly change the 
behavior investigated below. 

When is the assumption that <C 1 justihed? Recall that the characteristic radius 
Ps and radius for an NFW prohle are given in terms of the halo virial mass and 
concentration c: 

ta = CTs, ( 22 ) 




Ma = M(rA) = / p^Y^{r)dr = 

Jo 

Ps — ^cPcriti,^^- 


ln(l -I- c) 


l + C 


(23) 

(24) 


The density contrast 6c is in turn given by 


4 = 


A 

3"]^’ 


(25) 


where Kc = ln(l -|- c) — c/(l -f- c). The problem thus reduces to hnding an expression for A, 
the virial overdensity. In the spherical collapse model, this is given by A ~ for a 

hat universe [TOHTH] : A hence approaches the familiar value of 178 in the matter-dominated 
era. 


Inserting these expressions into (21) above yields an inequality for af in terms of c and 


Ma, along with the redshift of virialization z: 


af< -= (47r)-V3M-'/=* 

PsTs V 


f ^Pcrit ('^) 


- 2/3 


KcC 


-2 


(26) 


= 24.56 cm^/g x 


Ma 

IQi^M© 



Pcrit ( 2 :) 


Pcriti^Z 15) 



(27) 


In the second line we have inserted the typical halo parameters we will consider below: 
X = 15, Ma = IO^^Mq. 
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It remains to insert plausible values for the concentration c. Individual halos of mass 
Ma, formed at a hxed redshift z will have varying concentrations, but there should be some 
mass- and redshift-dependent median concentration, c{M/^,z). Prada et al. [H] used the 
Millennium |75U76j . Bolshoi m, and MultiDark [75] simulations to examine the shape of the 
c{M^, z) curve with varying mass and redshifts. They found that for each redshift considered 
(from z ~ 0 — 6) the concentration formed a U-shaped curve: it was minimized at a certain 
value of the mass, but increased steeply both above and below this mass. Furthermore, 
they found that both the minimum value of the concentration and the mass at which this 
minimum was realized decreased with increasing redshift. At the large redshifts we consider, 
the cluster-sized halos needed to form supermassive black holes are far more massive than 
the bottom of the U-shaped curve; accordingly, the htting formulae given in m predict 
that the concentration for these halos will be extremely large, of the order of c 10® for 
the halo parameters above. If this were true, the initial density prohles of these large, early 
halos would be extremely concentrated, so that their inner regions are extremely thick even 
for af > 10“®. In this case the simulations presented in this paper would not be reliable. 

We emphasize, however, that the htting formulae of m were devised using simulated 
halos only out to z ~ 6; they should not be trusted so far away from their domain of validity. 
Accordingly, we have consulted the high-redshift halo catalogs of the FIRE simulation [79] , 
which attempted to resolve an overdense region at high redshift. The catalogs use the Amiga 
Halo Finder [50] to measure c in the same way as dehned in 0. We are interested in the 
concentration parameters of the most massive halos formed at a given redshift. Perhaps 
unsurprisingly, we hnd that, even at z ~ 30, halo concentrations range from 2 to 11, similar 
to the values found at lower redshifts in the simulations consulted in na, rather than the 
much higher values predicted by naively applying the htting formulae. We do not attempt 
to construct the full c(Ma, z) curve at high redshifts on the basis of this limited data, but 
we do assume that realistic halos will take concentrations in this observed range. 

The upper bound on af for which < 1 ranges from 0.32 — 2.65 cm^/g as concentrations 
decreasing from 11 to 2 are inserted into (27). In the remainder of this paper we will typically 
set c = 9, which gives a bound of 0.425 cm^/g. In Section below we will hnd that this 
bound is of the same order of magnitude as the cross section needed to produce the desired 
high-redshift supermassive black holes using uSIDM. Accordingly, there is a surprisingly 
small region of parameter space where both the assumption of an initial NFW prohle is 
valid and the desired black holes are produced. We will discuss this further in Section 
1^ For now, we note only that the qualitative results of this paper should still hold even 
when our assumption of an initial NFW prohle is invalid. Outside of this range, we expect 
that gravothermal collapse should still occur—in fact, it should occur faster because core 
formation will have begun even before virialization—but the particular expressions given 
here will no longer be valid. 


C. Integration of the Equations 


Given the initial conditions, we can proceed to integrate the system of equations (12-17). 


We hrst move to a dimensionless form of the problem by choosing hducial mass and length 
scales {Mq, i?o}- Then the remaining dependent variables are given naturally in terms of 
these quantities, e.g. uq = GMq/Rq. Full expressions for all dependent variables in terms 
of Mq and Rq are given in Section 5 of [M]- The cross section per unit mass is now expressed 
dimensionlessly hj a = a/ao, ao = dvrRg/Mo. It is convenient to use the two quantities 
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FIG. 1: The dimensionless initial profiles for an NFW halo, with / = 0.01. The CDM and 
SIDM have the same velocity prohle, and their density prohles have the same shape but a 
normalization differing by //(I + /). As expected for SIDM, the initial luminosity at small 
radii is negative (shown as dashed on the plot), indicating that the cusp is being forced 
outward as a core begins to form. The glitch in the luminosity at r = 200 is a numerical 
artifact. 


already specihed in the NFW prohle, ’"s}; we therefore take 

{Mo = ^-kRIps, = g}- (28) 

Note that we have made a different choice of {Mq, i?o} than [M], since we consider a cuspy 
NFW prohle rather than a cored one and thus work with characteristic rather than central 
quantities. Finally, the timescale is set by the initial relaxation time at the characteristic 
radius, 

G,c(0) = l/^faps-UsCr), (29) 

so the independent variable can also be made dimensionless. Dimensionless quantities are 
written with tildes (e.g. p, t). The resulting initial prohles for an / = 0.01 halo are shown 
in Figure 
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We solve the problem by spatially discretizing into N concentric spherical shells, initially 
evenly logarithmically spaced in radius. At each timestep, we hrst apply the effects of heat 
conduction, which increases the energy within each shell, then adjust the prohle to maintain 
hydrostatic equilibrium. The heat conduction step is simple: we determine the luminosity 
prohle from the density and velocity dispersion using the dimensionless, discretized form 
of (15), then adjust the energy of each shell accordingly (using dU = Ldt for a hnite 
but small timestep). Timesteps are chosen so that the change of (dimensionless) specihc 
energy Ui = 3i>f/2 is not large: we require Aui/ui < e ^ 1 for each shell i, typically 
taking £ = 0.001. This means that as the gravothermal catastrophe approaches and core 
temperatures and densities become large, the size of timesteps will decrease dramatically: 
as expected, we cannot integrate through the collapse because the huid approximation itself 
breaks down there. 

To carry out the hydrostatic equilibrium step, we use the method of Lagrangian zones, 
in which the radius of each shell is adjusted while the mass it contains is left constant. 
The relaxation process, which involves long-range gravitational interactions rather than 
heat conduction via collisions, is entropy-preserving, so it preserves the adiabatic invariants 
Ai = piV^ ' for each shell i. After the heat conduction step, each shell is temporarily 
out of hydrostatic equilibrium, so that the equality (13) is violated by some amount Aj. 
The problem is to adjust the density, velocity dispersion, and radius of each shell i, such 
that hydrostatic equilibrium is again satisfied (Aj = 0 Vi) while preserving the adiabatic 
invariants. The assumption of adiabaticity, along with the use of Lagrangian zones to keep 
the mass of each shell hxed, hxes the density and velocity changes as functions of the set 
of changes of radii Ahj. Hence the requirement of hydrostatic equilibrium gives a system of 
differential equations for the changes of radii which, when linearized, is tridiagonal (since 
the thickness of each shell depends not only on its own central radius but that of its nearest 
neighbors). The resulting system is solved using a standard linear algebra library^. 


D. Results 


Unfortunately, the above procedure is still insufficient to integrate (12-17) in full gener¬ 


ality. The problem is that, because the SIDM and collisionless dark matter are separately in 
hydrostatic equilibrium, the method of Lagrangian zones will result in different sets of radii 
for the two species. But in order to perform subsequent timesteps, we need the total mass 
distribution at each radius for both types of DM. For computationally feasible numbers of 
shells {N ~ 400), interpolation is not accurate enough to preserve numerical stability and 
the distributions cannot be integrated all the way up to the point of gravothermal collapse. 

We can, however, consider the two limiting cases. (Luckily, these happen to be the cases 
we are interested in!) In the pure SIDM case / = 1, there is only one species and the problem 
does not arise. In the uSIDM case, / -C 1, we can ignore the gravitational backreaction 
of the uSIDM component on the collisionless DM and assume that it maintains an NFW 
prohle throughout, allowing the calculation of its mass distribution analytically at every 
point. We expect that the two cases should yield similar results, because the temporary 


violation of (13), the hydrostatic equilibrium condition, after each heat conduction timestep 


is overwhelmingly due to the increase on the LHS of the equation, from heat conduction. 


^ http://www.gnu.org/software/gsl, 
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(b) 


FIG. 2: (a) Evolution of SIDM density profiles, starting with an / = 0.01 NFW halo. Only 
the inner part of the halo is shown; the outer part still asymptotes to r~^ as in Figure [T] 
above for all halos. From top to bottom, prohles are at 0.0, 2.51, 4.79, 7.07, 9.35, and 
11.63 central relaxation times. Because G cc p~^, this corresponds to integrating for 
~ 1000 relaxation times in an / = 1 halo. However, comparison to the / = 1 results below 
suggests that the density prohle flattens in the same manner, just f~^ times slower: 
evidently the non-interacting dark matter has little influence on the central SIDM 
evolution, (b) Evolution of an / = 1 halo starting from NFW initial conditions. For 
clarity, only the inner portion of the density prohle is shown: the outer prohle has not yet 
changed signihcantly at this stage. From top to bottom, prohles are at 0.0, 9.73, 22.87, 
36.86, and 65.49 central relaxation times. As in the / = 0.01 case, the density prohle is 
battening as a core develops. 
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FIG. 3: Runaway collapse of an / = 1 SIDM halo with a = 0.088, starting from an initial 
NFW prohle. The inner prohle starts cuspy, rapidly shrinks to a self-similar prohle (as in 
[54] and Figure]^ above) with a p = 1 core, then slowly increases in density in a self-similar 
manner. After ~ 450 relaxation times, the core of the halo becomes optically thick, and 
self-similarity is broken: the core splits into a very dense inner core and an outer core 
which transitions between the two regions. Catastrophic collapse occurs as f/fr,c(0) 
approaches 455.65. 


rather than from interactions with the collisionless component, on the RHS of the equation. 
This is just the statement that the self-interaction is much larger than gravitational strength. 
We indeed hnd that this is the case, at least qualitatively. Consider Figure]^ which shows 
the early evolution of / = 0.01 and / = 1 halos with the same value of a. We see that 
behavior is indeed qualitatively the same: in both cases, a core begins to form as heat 
conduction dissolves the initial cusp. Note that the time scales are different: in the uSIDM 
case the relaxation time is increased by a factor of since the uSIDM density is a factor 
of / lower. So Figure suggests that uSIDM evolution is the same as the / = 1 case, just 
times slower. 

We will focus on the / = 1 case in the following, and then rescale our hnal results by f~^ 
as just described. Figures]^ andshow the entire evolution of an / = 1 halo with a = 0.088 
(chosen to allow comparison with |55]) from initial NFW prohle through to gravothermal 
collapse. First consider Figure]^ which shows the evolution of the density prohle. Although 
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Log Base 10 of Dimensionless Radius 


FIG. 4: Mass profile history of a cored SIDM halo with a = 0.088, starting from an initial 
NFW prohle. Once the core enters the optically thick regime, around t/tr,c(0) = 450, the 
inner core contains a constant total mass, around 2.5 — 3% of the characteristic mass Mq. 


the halo is initially in an NFW profile, the initial negative luminosity at small radii causes 
the cusp to empty out, driving evolution towards the cored, self-similar profile found by [H4] . 
as was already seen in Figure above. When the self-similar profile is reached after a few 
tens of relaxation times, the luminosity prohle becomes everywhere positive, and the core 
increases in density while its mass steadily shrinks. While the entire prohle is in the Imfp 
regime, evolution is self-similar, and the central density increases steadily. Inevitably, there 
comes a time, about 450 relaxation times after virialization, when the inner density increases 
enough that the most central regions enter the smfp regime, and the core bifurcates into a 
very dense outer core and an inner core which transitions between the two regions. 

Importantly, once the smfp regime has been reached, mass loss from the inner core is no 
longer efficient: the inner core has become so thick that evaporation is only possible from 
its boundary, not from the entire volume. This means that the mass in the inner core is 
essentially constant over the very short time (< 10G,c(0)) between breaking of self-similarity 
and catastrophic collapse. As mentioned in subsection 3 C above, the size of successive 
timesteps decreases rapidly as the gravothermal catastrophe approaches, so this short time 
takes very many (increasingly small) timesteps to integrate over, and the time of collapse 
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can be precisely given as 455.65 relaxation times after the start of integration. Becanse 
evaporation is inefficient after the loss of self-similarity, the mass in the inner core is still 
nonzero at the moment of collapse, nnlike in the globular cluster case, and a black hole will 
form. Figure 1^ shows that the inner core at collapse contains a mass of around 0.025Mo. 
Because the fluid approximation breaks down, we do not know that the entire inner core 
will collapse directly into a black hole, but, because it is optically thick, Bondi accretion [56] 
is extremely efficient. Hence, we expect that the black hole will rapidly grow to encompass 
the entire region regardless. 


4. SUPERMASSIVE BLACK HOLES FROM USIDM 

We have found that halos with a uSIDM component (and pure SIDM halos, on much 
longer timescales) grow black holes of mass Mbh = 0.025/Mo in a time 455.65tr,c(0). Given 
the considerations discussed above, uSIDM can help explain the existence of massive high- 
redshift quasars if the resulting black holes are large enough and form early enough that 
baryonic accretion can grow them to 10®Mq by 2 ; > 6. It remains to evaluate Mq and 
G,c(0) in terms of the halo parameters and use this requirement to place constraints on the 
uSIDM parameters {a, /}. 


A. Halo Parameters 


Instead of using the characteristic NFW parameters h is convenient to again 

parameterize a halo by its virial mass and concentration c. In dimensionless units, the 
mass contained within the zth shell is 

Mj = / pf'^df = / f~^ (1 -|- r'^df = - + ln(l -|- fj), (30) 

Jo Jo 1 + D 

but the virial radius ta = cvg, so 


Mbh _ Mbh _ 0.025/ 

Ma M{c) ln(l c) - c/(1 + c)' 


(31) 


This gives the desired expression for the seed black hole mass Mbh in term of the halo and 
uSIDM parameters. The denominator ranges from ~ 0.5 — 2 for realistic values of the halo 
concentration, so the BH mass is a few percent of the total uSIDM mass in the halo. 

Recall that the relaxation time is G,c(0) = l/{afpsi'sO'), i.e. the scattering time at the 
characteristic radius. The lower end of the interesting range for / = 1 SIDM is ~ 0.1 cm^/g, 
for which the relaxation time at the characteristic radius of a Milky-Way scale halo is 
approximately a Hubble time. To grow a black hole in galactic halos by z ~ 6, the relaxation 
time needs to be ~ 10 "^ times smaller to ensure ~ 500 relaxation times by the time the 
universe was a twentieth of its present age. This does not mean that a ~ 1000 /“^ cm^/g, 
though! Recall that ps = dcPcrit, where Sc is a function of c and the cosmology given below, 
and the critical density goes as (1 -1- z)^ in the matter-dominated era. Also oc ta oc 
{Ma/P critY^^ implies Ug oc Ma/ta oc Hence the mass and approximate redshift 

dependence of the relaxation time are 


( 32 ) 
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and we expect that af need not be that mnch larger than the interesting range for / = 1, 
i.e. we expect cr/ > 0.1 — 1 cni^/g. 

At this point, the reader might worry that this conclnsion combined with the observation 
of non-collapsed cores in the nearby {z ~ 0) nniverse rnles ont the existence of standard 
(/ ~ 1) SIDM. We emphasize, however, that large valnes of af mean that core collapse in 
times mnch smaller than the age of the nniverse is possible, bnt not, we expect, typical, it 
occnrs only in the rare halos which virialize at very high redshifts and remain nninterrnpted, 
i.e. do not experience major mergers, for long enongh to complete the gravothermal collapse 
process. See snbsection 5 A below for farther discnssion of this point. 

The exact expression for the halo relaxation time is 


tr,c(0) = 




afa 




(33) 


= 0.354 Myr x 


(- 


Ma 


VlO^^Mo 


- 1/3 


K, 


3/2 


/ Pcritjz) 


\ Pcriti^Z 15) 


- 7/6 




1 cm^/g 


-1 


(34) 

where Kc = ln(l+ c) —c/(l + c), 6c = {A/3)c^/Kc, and A, the virial overdensity, is 
for a flat nniverse, approximately 178 in the matter-dominated era. So the relaxation time, 
and hence the collapse time, is given in terms of the halo and nSIDM parameters. To match 
observations, we need some seed black holes to grow by a large enongh factor via Eddington 
accretion to reach Mbh ~ 10®Mq by 2 ; > 6; this leads to an ineqnality on a when the halo 
parameters and / are specihed. 


B. Explaining Observations 

Let ns spell ont the procednre more precisely. An observation of a particular high-redshift 
quasar at redshift Zohs yields a value for the luminosity, which corresponds to a supermassive 
black hole of mass Msmbh once the measured luminosity is identihed with the Eddington 
luminosity and a particular value for the radiative efficiency Cr is assumed. (We have already 
discussed potential issues with these assumptions in section]^ above; in the remainder of the 
paper, we will take the published observations at face value and assume their quoted SMBH 
masses, which take 6^=0.1 as input, are correct.) 

At the same time, the uSIDM framework developed in this paper tells us that NEW halos 
of viral mass Ma and concentration c virialized at redshift z form seed black holes of mass 
Mbh in a time 455.65tr^c(0), i.e. the seed black holes are formed at redshift Zcoii, where 


t{zcoii)-t{z)=455.65trA0), (35) 

and the time t{z) after the Big Bang corresponding to redshift z is given by the usual 
cosmology-dependent expression. 


t{z) = to 


da 

1 d ' 


(36) 


Equations (31) and (33) then give expressions for these quantities in terms of the halo 
properties {Ma, c, z} and the uSIDM parameters {a, /}. 
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There is still one parameter that must be specified: the fraction of SMBH mass which is 
due to accretion of baryons as opposed to the initial seed black hole. For simplicity, we will 
assume that the central black hole accretes continuously at the Eddington limit from the 
time of formation to the time at which it is observed. Of course more complicated growth 
histories are both possible and likely. Nevertheless, this simplifying assumption allows us to 
specify the fraction by instead giving N^, the number of e-folds of accretion at the Eddington 
limit. This finally allows us to compute the observable quantities: we have 


MsmBH = Mbh exp{Ne). 


(37) 

(38) 


To find acceptable values of a and / given the SMBH observables, we must specify (or 
marginalize over) the halo parameters and the baryonic contribution to the SMBH mass. 


The latter quantity directly sets (37) the redshift of seed black hole collapse, Zcoii, which 


yields the required collapse time and thus the required value of af via (33). Knowing the 


growth due to accretion of baryons also tells us (38) the required seed black hole mass Mj 


which specifies / via (31). 


'Sir, 


C. Examples 


As an example, consider again ULAS J1120-I-0641, with mass Msmbh ~ 2 x IO^Mq at 
Zobs = 7.085. To grow four orders of magnitude (A"e = In 10^) by Eddington-limited baryon 
accretion, for example, we must form a seed black hole with mass Mbh = 2 x IO^Mq by 
ZcoU = 12.9. With a halo of mass = IO^^Mq and concentration c = 9 formed at redshift 
X = 15, we find that tr,c(0) = 0.354 Myr x (1 cm^/g)/((j/). In order for 455.65 relaxation 
times to have passed in the 64.5 Myr between x = 15 and Zcoii = 12.9, we must have 
af = 2.50 cm^/g. From (31), we require / = 1.12 x 10“^ to get the correct seed mass, so 
cr = 2.23 X 10^ cm^/g = 3.97 x 10^ b/GeV. The large value of a is unsurprising: we chose 
to start with a halo much larger than the seed black hole we wanted to form, so / had to 
be small and a large in order to compensate. 

Alternatively, we could start with the same halo but produce the black hole entirely from 
uSIDM. The relaxation time is unchanged: ^, 0 ( 0 ) = 0.354 Myr x (1 cm^/g)/((j/). But 
now the black hole need not form until x = 7.085, 479 Myr after halo formation, so the 


required value of af is smaller, af = 0.336 cm^/g. Again applying (31) yields / = 0.112, 
a = 2.99 cm^/g = 5.36 b/GeV, coming much closer to the classic SIDM cross section. 

Of course, in the absence of direct measurements of the host halo of ULAS J1120-1-0641 
the problem is underdetermined. The point is that af takes reasonable values of 0{1) cm^/g, 
well within the regime described by the gravothermal fluid approximation starting from an 
initial NFW profile. 


5. DISCUSSION 

The examples in the previous section show that individual observations of high-redshift 
quasars can successfully be explained within the uSIDM paradigm. For uSIDM to be fruitful, 
however, we should ideally be able to find (or rule out) a consistent choice of these param¬ 
eters which successfully explains the cosmological abundance of high-redshift quasars. It is 
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unsurprising that some choice of cross section per unit mass a and fraction / can reproduce 
one particular observation, e.g. ULAS J1120+0641, but it is more suggestive if that choice 
can reproduce the entire observed number density of supermassive black holes as a function 
of mass and redshift. The minimal requirement for a viable uSlDM model is that it explain 
(or at least not conflict with) what has currently been observed. That means producing the 
correct abundance of ~ 10 ®Mq quasars at redshift 6 — 7, as has already been discussed, and 
ensuring that supermassive black holes are not overproduced in the nearby (lower-redshift) 
universe. Beyond that, one would like to make concrete predictions for the next generation 
of experiments, which should be sensitive to smaller masses and higher redshifts. 

This task is difficult for a number of reasons. The essential problem is that a number of 
nuisance parameters must be constrained or marginalized over in order to connect the uSlDM 
properties to the SMBH distribution (and then further to the quasar distribution). Even in 
the simplihed setup described above there were already the e-folds of baryonic accretion, N^, 
and the halo parameters Ma and c. In the cosmological context, these nuisance parameters 
are promoted to entire unknown functions that are currently only poorly constrained by 
observations and simulations. Even when constraints or functional forms are available, 
they are often trustworthy only in regimes far separated from the ones of interest to us 
here (for example, in the low-redshift universe, or in a lower mass range). We have already 
encountered this problem in Sectionabove, when considering the concentrations of massive 
NEW prohles at high redshifts. 

Nevertheless, in the remainder of this section we attempt to estimate the constraints that 
our existing knowledge places on the uSlDM parameter space. We hrst explain the source 
of our cosmological uncertainty and means by which it could be improved. Next we note 
a different source of tension within ACDM, independent of the existence of high-redshift 
SMBHs, that could be relieved by uSlDM. Finally, we present tentative maps of the uSlDM 
parameter space relevant to the resolution of these tensions. 


A. Cosmological Caveats 

Predicting the cosmological consequences of gravothermal collapse given a choice of the 
uSlDM parameters requires a unihed picture of the SIDM prohle at galaxy formation in 
terms of the halo mass and redshift, which will be easier given proper A-body simulations 
of halos containing uSlDM. There are several reasons why using the fluid approximation to 
simulate an isolated halo does not suffice. 

First, although the process of gravothermal collapse can be quite short on cosmological 
timescales, which is why it allows massive quasars to form faster than in the standard 
ACDM picture, we have seen that it is long in terms of halo time scales (several hundred 
characteristic relaxation times). It is therefore necessary for the halo to remain essentially 
undisturbed for this length of time in order for core collapse to occur and seed black holes to 
form.. The beginning and end of the collapse process—the elimination of the initial cusp and 
the catastrophic collapse itself after the core becomes optically thick—are driven entirely by 
dynamics in the innermost part of the prohle, so we might expect them to be insensitive to 
accretion or mergers in the outer halo. Figures and make clear, though, that these stages 
are very short compared to the length of the overall process. The vast majority of the time 
required for collapse involves the slow increase of density in the core as mass hows inward 
from the outer halo, which we expect to be sensitive to accretion or mergers. In other words, 
the halo must be isolated for several hundred relaxation times. Strong interactions with other 
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masses, such as major mergers, will disrupt the collapse process, essentially resetting the 
clock for seed black hole formation. Even accounting for more controlled accretion via minor 
mergers will technically necessitate the tracking of substructure within the collapsing halo, 
since it breaks the spherical symmetry required by the fluid approximation, although we 
expect it will not change our qualitative conclusions. Such tracking of substructure is only 
truly possible using A^-body simulations. 

More importantly, determining how often the collapse process is disrupted, and therefore 
predicting the spectrum of black hole masses as a function of redshift for particular values 
of a and /, in order to compare with existing and upcoming observations, requires detailed 
cosmological information. We need not only the halo mass function at very high redshifts 
(up to redshift 15 in the above example, and ideally out to at least 2 ; > 30 — 50) but also 
information on halo shape (the concentration parameter c(Ma, z) at the same high values 
of z, in the case that the halos form in NEW prohles) and, most importantly, detailed 
merger probabilities and histories as functions of mass and redshift. Even when analytical 
approximations to these quantities at z < 1 exist, it is unclear how conhdently they can be 
extrapolated to z ~ 50. Hence dedicated iV-body simulations are desirable. We will briefly 
note some additional interesting results, beyond the prediction of the history of the black 
hole mass function, which could be investigated given this cosmological information. 


B. The Too Big to Fail Problem 

This paper has noted that gravothermal collapse of uSIDM can produce seed black holes in 
the center of virialized halos. We have primarily been concerned with using this mechanism 
to explain the abundance of massive high-redshift quasars, but we now mention a few other 
areas where it could prove useful. We emphasize that these are logically independent of the 
quasar issue: we should not necessarily expect that the same choice of uSIDM parameters 
will be useful in both cases. 

First, it is intriguing that there exists a well-known (and relatively tight) relation between 
the properties of a host galaxy and the massive black hole it contains, the M-a relation 
[8Hj83] , which suggests some sort of causal mechanism connecting the central portions of the 
galaxy containing the black hole with the more distant regions where the velocity dispersion 
is measured. Gravothermal collapse naturally provides one such mechanism, and it would 
be suggestive if it produced the correct relation for some choice of the uSIDM parameters. 
At a minimum, it should not spoil the observed relationship in nearby galaxies; this has 
been used previously to constrain the cross section of / = 1 SIDM [HU ES] • 

More speculatively, the presence of central black holes in dwarf galaxies could resolve 
the “too big to fail” problem [HHl EZ] , in which the central densities of the brightest Milky 
Way satellites have much lower central densities than the most massive subhalos in ACDM 
simulations of Milky-Way sized galaxies. One way to resolve the problem is to invoke 
physics not present in the simulations to reduce the central densities (within ~ 1 kpc of the 
subhalo center) by a factor of order unity. If all of the dark matter is self-interacting with 
(T ~ 0.1 cm^/g, it naturally smooths out cusps to form cores, which could provide the needed 
reduction in density HD]. But the small fractions / 1 we consider in this paper cannot 

solve the problem in this manner; another method of removing substantial mass from the 
central ~ kpc is needed. 

Under some circumstances, it is possible that black holes could provide the needed re¬ 
duction in mass. Merging black hole binaries emit gravitational waves anistropically and 
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thus receive an impulsive kick, up to several hundred km/s. This energy can be distributed 
to the surrounding baryons and kick out a substantial portion of the central mass, forming 
a core gH EHl ES]. Such a scenario is only viable if the required binary black hole mergers 
are sufficiently common within dwarf galaxies or their progenitors. Although the standard 
cosmological model predicts the presence of black holes in the center of nearly all large halos, 
it is not clear that ACDM produces enough black holes within the smaller halos which are 
the progenitors of dwarf galaxies. Here we propose instead to use uSIDM to produce them. 

Solving the Too Big to Fail problem using black holes formed from uSIDM requires 
a particular sequence of events: hrst, small halos must remain isolated enough to form 
seed black holes; second, the probability of major mergers must become large enough that 
essentially all of the Milky Way satellites have binary black holes coalesce within them in 
order to reduce their central densities. During the epoch of matter domination, we see that 
the black hole formation time for a halo of hxed mass goes as (1 + (32), while we 

expect the merger timescale to be set roughly by the Hubble time, H -\z) ~ (l + z)-3/^ 
So halos of a given mass that form before some critical redshift will indeed grow black holes 
before they merge. In the next subsection we consider the parameter space where black hole 
seeds are ubiquitously formed in the progenitors of today’s dwarf galaxies. 


C. Parameter Space 


1. High-Redshift Quasars 


In subsection |4 Cj above, we presented two possible routes to produce a supermassive black 
hole matching observations. Here we move from specihc examples to a discussion of the entire 
parameter space relevant to the production of high-redshift quasars like ULAS J1120+0641. 
Recall that we have six input parameters: {M^,c,z,Ne,o', f}, respectively the halo mass, 
concentration, redshift of virialization, e-folds of Eddington-limited accretion after collapse, 
and uSIDM cross section per unit mass and fraction. We specify the halo properties as 
above: Ma = lO^^M©, c = 9, z = 15. We then use ( |^ and the redshift at which a quasar 
is observed, in this case Zobs = 7.085, to eliminate Ne, leaving a two-dimensional parameter 
space for production of black holes by this time. Finally, the requirement that the mass 
of the quasar match observations, Msmbh ~ 2 x IO^Mq, combined with the assumption 
of continuous Eddington-limited growth since black hole formation, reduces the parameter 
space to one dimension, a curve cr(/). 

We present the parameter space in Figure The one-dimensional curve cr(/), where 
continual Eddington-limited accretion since black hole formation results in a supermassive 
black hole with Msmbh = 2 x IO^Mq at Zobs = 7.085, is the solid black line. Because 
the baryonic accretion history after seed black hole formation is uncertain, as discussed in 
Section above, we also indicate with the shaded regions the entire portion of the full a~f 
plane in which black holes of any size smaller than Msmbh are produced by Zobs- 

There are several constraints on this reduced parameter space. First is the simple require¬ 
ment that gravothermal collapse indeed occurs before Zobs = 7.085. We have already seen 


in subsection 4C above that this constrains af > 0.336 cwr/g. Second is the requirement 


that the black hole produced by gravothermal collapse must not be larger than the observed 
mass of ULAS J1120+0641. Combined with additional assumptions about baryonic accre¬ 
tion, this excludes the entire region above the black curve in Figure [5j Even without the 
assumptions, this still constrains the black hole mass via equation (31), and therefore the 
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FIG. 5: uSIDM parameter space for production of massive high-redshift quasars. We have 
used the numbers considered in the example above: Msmbh ~ 2 x Zobs = 7.085, 

Ma = c = 9, X = 15. The solid line plots values of a and / that result in an 

SMBH of the desired size at the time of observation, assuming continuous Eddington 
accretion from the time the core collapses and the seed black hole is formed. The green 
dotted vertical line marks the largest allowed value of /. To its right, collapsed black holes 
are already larger than Msmbh- To its left, collapsed black holes form smaller than 
Msmbh, but can grow larger by accreting baryons. The points on the bine dashed line all 
result in collapse precisely at the redshift of observation; below this line, a black hole has 
not yet formed by Zobs- Points on and above the red dashed line result in a halo that is 
already optically thick at the time of virialization, i.e. optically thick at the characteristic 
radius for the initial NFW prohle. As discussed in the text, the methods used in this paper 
are not directly applicable here, bnt we still expect gravothermal collapse. Numerical 
values for all of these bounding lines are given in the text. 


uSIDM fraction /, provided that a black hole is actually produced. Again, the resulting 
constraint was calculated in snbsection |4 C[ / < 0.112. Larger valnes of / would produce 
black holes which contained too large a portion of the mass of the entire halo. 

Finally, recall that onr expressions for the collapse time and resnlting black hole mass 
are based on simulations. As discussed in snbsection 3 B above, the simnlations assnme 
the nSIDM is initially in NFW prohle, which is only valid if nSIDM interactions were slow 
compared to the timescale of halo formation, i.e. when the halo is initially optically thin. 
This places a constraint on af as a fnnction of the halo parameters, given by eqnation (27). 








23 


For our chosen values this gives af < 0.425 cm^/g. 


Because af directly sets the collapse time via (33, 35), the upper bound on af is also 
a lower bound on the time of formation of a black hole from an initially optically thin 
uSIDM halo: in this case, we must have Zcoii < 7.90. In turn, this places an upper bound 
on the number of e-folds of growth from baryons that can occur before Zobs = 7.085, via 
We hnd W < 2.24, i.e. black holes formed from optically thin uSIDM halos have time 
to grow less than an order of magnitude from baryons. In particular, we cannot trust the 
precise results of our simulations for the example we considered in subsection 4 C above, with 
Ne = In 10^. Note, however, that the upper bound on af, and thus on Zcoii, is independent 
of Zobs- it depends only on z, the redshift of halo formation. Most high-redshift quasars are 
seen near Zobs ~ 6: ULAS J1120-1-0641 is an outlier. Black holes at redshift 6 have had time 
for another 180 Myr ~ 4tsai of baryonic growth, so they can grow by up to a factor of ~ 500 
from baryons. 

We have seen that there is an extremely narrow range, 0.336 cm^/g < af < 0.425 cm^/g, 
in which the uSIDM halo considered here is optically thin at virialization but nevertheless 
rapidly collapses to form a black hole. How can we explain the closeness of these two 
bounds? In general, they are not independent. The upper bound requires that the initial 
halo be optically thin, i.e. that the scattering cross section be less than the “characteristic 
cross section” of the halo, l/(ps^s)- But the lower bound requires that collapse not take too 
long, i.e. that the scattering time at the characteristic radius, l/iafpsi's) is small compared 
to a Hubble time. These bounds can be simultaneously satished when ~ f's/us. But, 
ignoring concentration dependence and numerical factors, Ug ~ ^/GM^JVs ~ ~ 

Hts, so rs/vs ~ as desired. That is, both bounds exhibit the same mass dependence, 
and their redshift dependence is identical when z and Zcoii are similar, as can be verihed from 
(26, 33, 35). For our particular choice of c = 9, the numerical factors are nearly canceled by 
the concentration dependence, so the bounds are especially close. 

We emphasize again, however, that the upper bound on af (the red dashed curve in 
Figure is not a true physical exclusion of the uSIDM parameter space above it. It 
merely signals that the fluid approximation used in this paper is no longer valid outside 
this space. As discussed extensively in subsection 3^ above, we expect that prohles in 
which the uSIDM starts optically thick should in fact undergo collapse even faster. The 
requirement of an optically thin initial prohle would only be physical if starting otherwise 
led to fragmentation, turbulence, or some other mechanism by which the core was destroyed 
or core collapse avoided. 

Figure [^presents the uSIDM parameter space for a particular choice of halo parameters 
{Ma,c, 2;}. We briefly consider how constraints on the parameter space are changed when 
these parameters are altered. First consider the halo mass M^. The collapse time (33) scales 
as ' , so smaller values of the halo mass require values of af to form black holes in the 
same time. At the same time, the value of af required for the halo to start initially optically 
thin (26) has the same scaling with halo mass. So decreasing Ma will shift both bounds 
to higher values of af (up and to the right on the a-/ plane), but it will not qualitatively 
change the shape of the allowed parameter space. This shift accounts for the main difference 
between the high-redshift quasar parameter space and the dwarf satellite parameter space we 
will consider next. We note, however, that at relatively low redshifts there is a well-known 
black hole-bulge relation |8ll [90l |9T] , Msmbh ~ ^0~^Mbuige- If this relation persists at high 
redshifts, we should not depart too far from Ma ~ IO^^Mq to explain Msmbh ~ IO^Mq. 

Next consider the concentration parameter c. Again consulting (33), we see that the 
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collapse time depends strongly on concentration, scaling roughly as The collapse time 

depends more strongly on the concentration than does the optically thin condition, so for 
small enough values of c it will be impossible to form black holes before a given redshift 
starting from an optically thin halo. When the other halo parameters and Zobs are kept 
hxed, we hnd that this critical value of c is 7.4. Conversely, by going to larger and larger 
values of c we can form black holes by any desired time at smaller and smaller values of af. 
However, extremely high values of the concentration parameter correspond (unsurprisingly) 
to extremely concentrated halos, with <C r^,. It is not clear that such halos are actually 
produced in ACDM. 

Finally, consider the redshift of halo formation 2 ;. Once more consulting (33), we see 
that we can take the collapse time to zero by increasing z. Heuristically, this is because 
the critical density, and thus characteristic density, increases with increasing redshift, so a 
just-virialized halo is closer to the densities needed to start the catastrophic collapse process. 
However, producing large virialized halos at higher and higher redshifts becomes increasingly 
unphysical given the bottom-up structure formation mechanism in ACDM. Decreasing 2 ; has 
the opposite effect: halos of a given size become more common, but larger values of af are 
required to produce a black hole by a given Zobs- Like the case of small concentration 
parameter, for small enough ^ it is impossible to form black holes starting from optically 
thin halos before a given time. In this case the bound on the redshift of formation for the 
halo considered here is z > 13.53. 


2. Dwarf Galaxies 


Recall from the previous subsection that one resolution to the Too Big to Fail problem is 
the formation of cores in dwarf galaxies if matter is ejected during binary black hole mergers. 
Our goal here is to specify the parameter space in which uSIDM produces black holes in the 
progenitors of dwarf galaxies before the epoch in which binary mergers are common. As in 
the case of high-redshift quasars above, we start by specifying a set of typical values for the 
halo parameters {M^,c,z}. Ref. [87] compared the Milky Way dwarf galaxies to subhalos 
around similarly-sized galaxies in the Aquarius simulations [92| to derive probable values for 
the virial mass and maximum central velocity Umax of each halo at the time of its infall 
into the main Milky Way halo. 

Recall that the maximum velocity of an NFW prohle is 


^'max = 0.465 Ua, 



with Kc = ln(l -|- c) — c/(l -f c), at radius 

r^ax = 2.163rs, 


(39) 


(40) 


as can easily be verihed numerically using the dehnition of the NFW density prohle (19) 
and V = ^yGM{r)/r. Rearranging gives an expression for K^jc in terms of Ma, Umax, and 
Zinfaii- In particular, since r/^ ps ^ , we have Kdc ~ pcritdY^^- But Kcjc has 

a maximum value of 0.216, so above some value of Zinfaii, there is no possible NFW prohle 
with the given values of Ma and Umax- For the derived values for the Milky Way dwarfs, 
we hnd in general that Zinfaii ^ 6. As the infall redshift moves lower for each particular 
dwarf, the concentration increases from a minimal value c = 2.16. (Actually, K^jc attains 
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FIG. 6: uSIDM parameter space for production of black holes in dwarf galaxies. The 
parameters used are Msmbh = IO^Mq, Zobs = 4.5, Ma = 10®Mq, c = 9, z = 15. The solid 
line plots values of a and / that result in an SMBH of the desired size at the time of 
observation, assuming continuous Eddington accretion from the time the core collapses 
and the seed black hole is formed. The green dotted vertical line marks the largest allowed 
value of /. To its right, collapsed black holes are already larger than Msmbh- To its left, 
collapsed black holes form smaller than Msmbh: but can grow larger by accreting baryons. 
The points on the blue dashed line all result in collapse precisely at the redshift of 
observation; below this line, a black hole has not yet formed by Zobs- Points on and above 
the red dashed line result in a halo that is already optically thick at the time of 
virialization, i.e. optically thick at the characteristic radius for the initial NFW prohle. As 
discussed in the text, the methods used in this paper are not directly applicable here, but 
we still expect gravothermal collapse. Numerical values for all of these bounding lines are 
given in the text. 


its maximum at c = 2.16 and approaches zero both as c —)• 0 and c —>■ oo, but we neglect the 
former branch, with c < 2.16, as unphysical.) In particular, we choose Zinfaii ~ 4.5, which 
results in c ~ 9 for a typical dwarf, the same as considered above. 

A typical Milky Way dwarf in [H7| has Ma = 2 x 10®Mq at the time of infall. If Too 
Big to Fail is to be explained by means of binary black hole mergers, a typical dwarf should 
have undergone a major merger, so that a binary black hole merger occurs in the hrst place. 
We will therefore take Ma = 10®Mq, c = 9, Zobs = 4.5 as our typical parameters. 

Figure presents the parameter space for our typical dwarf halo. We have taken the 
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black hole mass to be Msmbh = lO^M©, in accordance with the black hole-bulge relation, 
and again assumed that the redshift of formation of the progenitor halo is ^ = 15. The 
various bounds in the hgure are attained in the same manner as they were for the quasar 
bounds shown in Figure]^ so we simply quote them here and refer to the discussion above 
for details of their calculation. The lower bound on af, which comes from requiring collapse 
before Zobs = 4.5, is af > 3.26 cm^/g. The upper bound on /, which is calculated using (31) 
and scales linearly with Msmbh, is / < 0.056. 

The upper bound on af, set by the requirement of an optically thin initial prohle, is 
o'f < 9.16 cm^/g. This corresponds to an upper bound on the redshift of collapse, Zcoii < 
7.90. (As discussed above, the upper and lower bounds have the same mass dependence, 
and we are considering the same values of 2 ; and c as we did for high-redshift quasars, so 
we recover the same bound on the collapse time.) Once again, this gives an upper bound 
on the number of e-folds of growth from baryons, via (37). Because we have a much longer 
time for growth after black hole formation than we did in the high-redshift quasar case, it is 
much looser: Ne < 15.2. Once again, the allowed range on af is a factor of only a few. But 
the signihcantly looser bound on Ng means that black holes can grow by a factor of nearly 
4 X 10®. So / can be decreased by over six orders of magnitude from its maximal value, 
and a increased by a corresponding amount, while still maintaining an optically thin initial 
prohle and allowing reasonably large black holes to form. This explains the much larger 
range in a and / seen on Figure [^compared to Figure]^ 


3. Both Simultaneously? 


We have just seen that the minimum value of af needed to produce massive high-redshift 
quasars is about an order of magnitude lower than those needed to produce black holes 
in dwarf galaxies before major mergers. We can understand this qualitatively from the 
expression for the halo relaxation time, (33): it scales as . The black holes in dwarf 


galaxies have about twice as long to form, until Zobs = 4.5 instead of Zobs = 7.085, so af is 
scaled by a factor of /2 10. 

This scaling of af implies that the uSIDM parameters which produce black holes in 
dwarfs are a strict subset of those which produce high-redshift quasars. It is then easy to 
choose values which solve both: one simply takes af > 3.26 cm^/g and chooses compatible 
values of a and / to taste. Because af is signihcantly larger than the minimum value 
needed to produce high-redshift quasars, the uSIDM halos which produce them will start 
initially optically thick, above the (red dashed) upper bound on af shown in Figure In 
this optically thick regime, the expressions for the gravothermal collapse time and black hole 
seed mass derived from the simulations of Section should be taken as limits: we expect 
that gravothermal collapse should occur a slightly shorter time after halo formation and 
result in slightly larger seed black holes. 

As an example, consider the {a, /} values that fall in the one-dimensional parameter 
space discussed at the beginning of this subsection, where continuous Eddington accretion 
from the time of black hole collapse until Zobs just produces a supermassive black hole with 
the observed value of Msmbh = 2 x 10®Mq (the solid black line in Figure]^. If af = 3.26, 
the smallest possible value needed to also produce black holes in dwarfs by redshift 4.5, we 


would hnd using (33, 37) that Zcoii = 13.3. (Strictly speaking (37) is not valid in the context 


of an optically thick initial halo, since it uses an expression for the collapse time derived 
from an initially optically thin NFW prohle. We are simply using it here for the sake of 
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illustration.) In this case there is time for 9.54 e-folds of baryonic accretion before z = 7.085, 
and the initial black hole has mass 7 x IO^Mq. This determines the USDIM parameters for 
this example as cr = 8.14 x 10® cm^/g, / = 4.01 x 10“®. 

As the beginning of this section emphasized, it is largely beyond this scope of this paper 
to describe a fully consistent uSIDM cosmology. Nevertheless, this subsection suggests that 
the USDIM paradigm is flexible enough to resolve both of the potential tensions within 
ACDM discussed here. It is possible that a single species of ultra-strongly self-interacting 
dark matter could in fact resolve both tensions simultaneously. In investigating this question 
further, it will be important to move beyond the simplifying assumptions employed herein, 
especially the stipulations of an initial optically thin prohle and a cosmologically isolated 
prohle. 


6. CONCLUSION 

In this paper, we considered a minimal extension of the SIDM parameter space, in which 
a self-interacting component comprises only a fraction of the dark matter. For / < 0.1, 
this evades all prior constraints on SIDM models. We highlighted the uSIDM regime, where 
the SIDM component is subdominant but ultra-strongly self-interacting, with / ^ 1 and 
(j ^ 1 cm^/g. In the setup considered here, the presence of uSIDM leads to the production 
of black holes with a mass of around 2% of the total uSIDM mass in the halo at very 
early times. In particular, such black holes can act as seeds for baryon accretion starting 
soon after halo formation, alleviating potential difficulties with accommodating massive 
quasars at high redshifts within the standard ACDM cosmology. If black holes are formed 
ubiquitously in dwarf halos before they undergo mergers, they may also resolve the Too Big 
to Fail problem by ejecting matter from cores during black hole mergers. More detailed 
cosmological simulations are needed to confirm the conclusions of this paper and suggest 
other potential observational consequences of uSIDM. 

Setting aside the detailed predictions, this paper has demonstrated that multi-component 
dark matter can have strong effects on small scales while still evading existing constraints. In 
the toy model discussed here, the strong effect was the result of the gravothermal catastrophe. 
Gravothermal collapse of a strongly-interacting dark matter component is a novel mecha¬ 
nism for production of seed black holes, potentially one with many implications. Given its 
appearance in the simple extension of ACDM considered here, it is plausible that gravother¬ 
mal collapse and its observational consequences, such as seed black hole formation, are 
generic features of more detailed models. It is important to consider, and then observe or 
constrain, this and other observational consequences that are qualitatively different from 
the predictions of the standard cosmological model. 

Our discussion has been purely phenomenological, so it is reassuring to note the existence 
of a class of hidden-sector models [93] which naturally produce a subdominant strongly- 
interacting dark matter component, with self-interaction cross-sections ranging as high as 
a ~ 10^^ cm^/g. Very interestingly, some models give both a dominant component with 
(T ~ 0.1 — 1 cm^/g, as needed to alleviate discrepancies between ACDM and observations, 
and a uSIDM component with a ~ 10® — 10^ cm^/g, which could produce seed black holes 
via the mechanism described in this paper. 
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